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ABSTRACT 

Equations are presented for calculating power law exponents from wind 
speed and surface roughness data. Results are evaluated by comparison with 
wind profile date measured at a variety of sites. 

INTRODUCTION 

The power law equation is a simple yet useful model of the vertical 
wind profile which was first proposed by Heilman (1916), according to Simiu 
and Scanlan (1978) . The general from of this equation is 

V2 = ViCz2/2i)®^ CD 

in which Vj and V 2 are simultaneous steady wind speeds over level ter- 
rain at elevations zi and 22 , respectively. The exponent OL is deter- 
T^ined experimentally. For example, early work by Von Karman (1921) showed 
that under certain conditions DC is equal to 1/7, indicating a correspon- 
dence between wind profiles and fluid flow over flat plates (Schlicting, 
1968). In the general case, however, 0(. is a highly variable quantity. 
Sisterson and Frenzen (1978) measured wind profile exponents which changed 
from 1/7 during the day tc 1/2 at night, over the same terrain. Golding 
(1955) describes OC as an exponent which varies with height, time of day, 
.season of the year, nature of the terrain, wind speed, and temperatiire . 

Most investigators agree that a constant value of is an oversimplifi- 
cation, and that OC must be treated as a statistical parameter. 

In spite of the variable nature of the exponent ft , Equation (1) is 
simple enough in form to permit adding some of the factors listed by Golding 
without losing the engineering convenience of the power law. Two of these 
factors will be considered here: (1) the nature of the terrain, in terms of 

its surface roughness, and (2) the wind speed. 

Up until now, the effects of wind speed and surface roughness have only 
been considered separately. As shown in Figure 1, several functional rela- 
tionships have been proposed for the variation of the exponent tt with wind 
speed. These functions are the constant 1/7 law, a step function (Fales, 
1967) , a linear function for wind speeds exceeding 27 meters per second 
(ASCE, 1961), » power function (Fichtl § Smith, 1977), and a logarithmic 
function (Justu:: g Mikhail, 1976). All of the variable functions show a 
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decrease in the exponent with increasing wind speed, but none includes the 
effect of surface roughness. Frost et al (1978) and Justus (1978) give 
equivalent values of the exponent and the well-known surface roughness 
length Zq, but without specifying the wind speed. 

Thus, information on the separate effects of wind speed and surface 
roughness is available, but a wind profile model which combines these two 
effects is lacking. The objective of this paper, then, is to present equa- 
tions relating the mean value of DC to both the surface roughness and the 
steady wind speed. Other factors are assumed to be responsible for statis- 
tical variations about this mean value. 

The equations presented in this paper are exploratory in nature and 
should be reevaluated with more data, when such data are available. 

SUMMARY OF EQUATIONS 

The Justus-Mikhail equation shown graphically in Figure 1 was selected 
as a starting point for this study because of its relative simplicity and 
the fact that it has been shown to be compatible with wind distribution 
data. This equation for the wind profile exponent DC in Equation (1) can 
be written as follows: 

DL = 0.37 - 0,203 log Vr (2) 

in which Vr is the wind speed at a reference evaluation of 10 meters. 

Equation (2) can be generalized to include a surface roughness para- 
.•neter and statistical variations. Assum^ OC is a random variable which is 
normally distributed about a mean value flC with standard deviation U' . 
Assume further that both_^ and qr are functions of the steady (non- 
gusting) wind speed and jZ. is_also dependent on the surface roughness, 
ihen, a proposed equation for , based on Equation (2) is 


1 

and 

1 - log Vi/log Vjj 
1 - Q(.o log Czi/tj.)/log Vh 

C3a) 

tX-o 

/■ / ^0.2 
* (Zo/Zj.) 

(3b) 

in which 




surface roughness exponent 


Vh 

homogeneous wind speed ( 5L ® 0) , m/s 


Vi 

steady wind speed at elevation Zj, m/s 


Zx 

reference elevation, 10m 



surface roughness length, m 
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Equation (3b) is an empirical relationship which will be discussed 
later. The standard deviation (T^(Vi) oust be determined from site data. 

Equations (3) are s*'Own graphically in Figures 2 and 3. The terrain 
descriptions in Figure 3 follow Frost et al (1978) and Justus (1978) . 

Weibull distribution parameters C and k and the wind profile para- 
meters and Vh are related, as shouu by the following equations: 


P(Vi>«V) 


exp 




(v/cp^'^j 


C4) 


in which 


P(Vi>V) 


Cl 


probability that the steady wind speed at elevation 
2 j will exceed V 

Weibull scale factor at elevation m/s 


kl 


Weibull shape factor at elevation Zj 


Weibull distribution parameters at an elevation Z 2 can be calculated 
from those at elevation z\ by means of the following equations: 



Cl (^2/^1^ 


Oi 


c.l 


a, 


1 - log Cj/log Vj^ 

® 1 - Uj/tp/log Vj, 


(Sal 

(5b) 


1 - C(o log (Zj/zp/log Vj^ 

1 - log (z 2 /Zy)/log Vj^ 


(5c) 


The Justus-::ikhail equation, (2), is a special case of Equation (3a) 
with SLo equal to 0.37 and equal to 67 meters per second. Eqtiations 
equivalent to Equations (3a) ar.d (5) with these specific values of bio and 
Vh were also used by Justus et al (1976) to calculate reference wind speeds. 

CALCULATED AND OBSERVED WIND PROFILES 
Mean Wind Speed Profiles 

Wind profile data (WSSI, 1976 to 1978) from wind turbine sites or poten- 
tial sites selected by the Department of Energy were used for a preliminary 
evaluation of Equations (3). Table I lists calculated roughness exponents 
for five such sites and also for the Justus-Mikhail reference site. 

A? homogeneous wind speed of 67 meters per second was assumed for all sites 
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in this prelieinary evaluation. Figure 4 illustrates the procedure for 

calculating , using data for Clayton, MM. Exponents calculated from 

monthly average wind speeds at several elevations are plotted versus the 

logarithm of the reference elevation wind speed, Vp. A straight line 

through the centroid of, these data points and the homogeneous point 

(Vr « Vh « 67 m/s and tC ” 0) defines . at a wind speed of one meter per 

second. 

The roughness exponents in Table I vary from 0.1 to almost 0.5, in 
qualitative agreement with the site characteristics. Higher values 
occur at rougher sites and lower values at smoother sites. The Justus- 
Mikhail constant of 0.37 appears to be somewhat too large to be used as an 
average roughness exponent. However, more profile data would have to be 
analyzed before a better reference value can be established. Based on these 
preliminary calculations, an average of about 0.30 is indicated. 

Equation 3(b) is an es^Jirical relationship between CCp and Zq, which 
is the conventional measure of surface roughness. As shown in Figure 5, the 
exponent of 0.2 was obtained by curve-fitting equivalent values of OC and 
zo from the literature (Frost, et al, 1978, and Justus, 1978). Calculated 
values of Table I then established the coefficient as (ZiO”®*^, 

Peak Wind Speed Profiles 

As shown in Figure 2, the proposed profile model is based on the assump- 
tion that wind profiles become uniform or "homogeneous" at high wind speeds, 
far all values of surface roughness. This concept represents an extrapola- 
tion of the Justus-Mikhail equation as originally proposed, and additional 
varification is required. To do this, profile exponents were calculated for 
monthly peak winds measured at 11 DOE sites (WSSI, 1976 to 1978), to deter- 
mine if there was a correlation between high winds and low profile exponents. 

Typical results of the peak wind analysis are shown in Figure 4, for 
the Clayton MM, site. Values of OC calculated from monthly peak winds 
which averaged 23 meters per second were found to be significantly lower 
than exponents calculated from monthly average winds of about 6 meters per 
£ccond. The assumed homogeneous wind speed of 67 meters per second is con- 
sistent with the Clayton peak wind data. 

Table II summarizes the analysis of peak wind profile data. This table 
lists mean peak wind speeds, mean profile exponents, and standard deviations, 
for each of 11 DOE sites. A co^i^osite average is also given which indicates 
that mean profile exponents less than 0.07 can be expected at wind speeds 
higher than 22 meters per secp’^d. To predict extreme wind loads, the "mean 
plus 33*" value of (jC, would be used, which is 0.27 for the composite average 
site. This compares closely with the ASCE recommended value of 0.30 shown in 
Figure 1, for a wind speed of 22 meters per second. Thus, equations (3a) and 
(3b) are representative of the mean values of the wind profile exponent, and 
should not be used in extreme winds analysis. More conservative values should 
be used. 
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Mean profile exponents in Table II are shown graphically in Figure 6. 
Also shown in this figure are mean profile exponents calculated from monthly 
average wind speeds measured at five DOE sites (see Table I). A band repre- 
senting predicted values of R for flat, open terrain is also shown for 
comparison. At low wind speeds surface roughness is an important considera- 
tion in determining Exponents for the two rough sites (Plum Brook and 
Block Island) fall above the band, and that for e smooth hillside site 
(Culebra) falls well below. At high wind speeds surface roughness is no 
longer a significant variable and mean exponents for all sites are low, 

CONCLUDING REMARKS 

Equations have been presented and evaluated for a wind profile model 
which incorporates both roughness and wind speed effects while retaining the 
basic simplicity of the Heilman power law. Moreover, these equations recog- 
nize the statistical nature of wind profiles and are compatible with exist- 
ing analytical models and recent wind profile data. 


Predictions of energy output based on the proposed profile equa- 
tions are 10 percent to 20 percent higher than those made with the 1/7 power 
law. In addition, correlation between calculated and observed blade loads 
is significantly better at higher wind speeds when the proposed wind profile 
model is used instead of a constant power model. 

It is recommended that statistical analysis of site wind data include a 
bivariate distribution of profile exponent versus wind speed. This would 
permit more accurate calculation of the parameters in the proposed model and 
permit estimates of standard deviation as a function of wind speed. Statis- 
tical information of this type would be extremely useful to wind turbine 
desigiers for predicting wind load spectra, fatigue life of components, and 
energy output. 
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Table I. - Surface roughness exponents for five candidate wind 
turbine sites. 


Site 


Surface roughness ® 
Type of Terrain exponent, 


Plum Brook, OH flat; high woods 0.48 

Block Island, RI hilly; medium woods 0,44 

Reference (Justus and relatively flat 0.37 

Mikhail, 1976) 

Clayton, NM flat; open 0.30 

Russell, KS flat; open 0.30 

Culebra, PR shoreline hillside; 0.10 

bushes 


Ot = - log Vy/log Vh) , with Vh = 67 m/s 


Table II. - Wind profile exponents for monthly peak winds at 
candidate wind turbine sites. 


Site 

No. of 
readings 

Mean peak 
speed at 
10m, m/s 

Wind profile 
Me^, 

R 

exponent , OC 
Std. dev., 

(T 

Amarillo, TX 

26 

22.3 

0.053 

0.045 

Block Island, RI 

30 

19.5 

.114 

.052 

Clayton, NM 

24 

22.8 

.054 

.044 

Cold Bay, AK 

9 

24.1 

.042 

.044 

Culebra, PR 

30 

17.4 

.010 

.078 

Huron, SD 

28 

20.6 

.094 

.056 

Kingsley Dam, NA 

29 

24.6 

.061 

.073 

Medicine Bow, WY 

5 

29.0 

.077 

.027 

Point Arena, CA 

23 

20.1 

.072 

.144 

Russell , KS 

29 

23.4 

.059 

.040 

San Gorgonio, CA 

33 

23.5 

.088 

.053 

Composite Site 

266 

22.2 

0.067 

0.068 
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Figure 2 - P^OP^ equation for wind profile exponent 

in terms of surface roughness and wind speed. 
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Figure 4 - Comparison of observed and calculated wind profile 
exponents for Clayton, New Mexico site. 
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Figure 5 - Curve-fit of data relating wind profile exponent tC 
to surface roughness length. 
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Figure 6 - Comparison of observed and calculated wind profile 
exponents for various sites. 


